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We consider the exclusive decay B -^ J/ipK using the QCD-improved factorization 

^ I method in the heavy quark Hmit. It is shown that the decay amphtude is factorizable 

^r ■ in this hmit and nonfactorizable contributions are calculable from first principles 

in perturbation theory. Also the spectator contributions at order a^ are finite and 



0^ I suppressed in the heavy quark limit. We present the result at next-to- leading order 

^? ■ in strong interaction, and leading order in l/m-b in the heavy quark limit. 

O; 

Q^l Exclusive nonleptonic decays of B mesons have received a lot of attention since 

they are observed in experiments at CLEO, BABAR and BELLE [1]. They 
offer the opportunity to test the unitarity of the Cabibbo-Kobayashi-Maskawa 
(CKM) triangle, the CP violation and to observe new physics effects. However, 
rS I quantitative understanding of nonleptonic decays of B mesons is difficult since 

c^ I there are final-state interactions which are intrinsically nonperturbative in 

nature. But when a B meson decays into two light mesons, there is large 
momentum transfer to final-state mesons. In this case, we can systematically 
calculate the decay rates using perturbation theory in the heavy quark limit. 

Beneke et al. [2] considered general two-body nonleptonic decays of B mesons 
extensively including a light-light meson system as well as a heavy-light system 
in the final state. The general idea is that in the limit rrih ^ Aqcd, the hadronic 
matrix elements can be schematically represented as 



{M^M2\0\B) = (Mi|ji|5)(M2|j2|0) 

X l + ^r„a^ + 0(AQCDM) 
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where Mi, M2 are final-state mesons and O is a local current-current operator 
in the weak effective Hamiltonian. If we neglect radiative corrections in a^ 
and power corrections in Aqcd, we get the factorized result with a form factor 
times a decay constant. At higher order in a^, this simple factorization is 
broken, but the corrections can be calculated systematically in terms of short- 
distance Wilson coefficients and meson light-cone distribution amplitudes. We 
call this the QCD-improved factorization. 

When we use light-cone meson wave functions for exclusive decays, the tran- 
sition amplitude of an operator Oj in the weak effective Hamiltonian is given 

by 



1 
{J/^K\0,\B)=J2Ff-''{ml) J dxTl^{x)M^) 

J 

1 

+ d^dxduT-\^,x,u)(j)Bi04''4>{x)4'Kiu), (2) 





where F^~*^ {m'^) are the form factors for B ^ K, and 0m (a;) is the light- 
cone wave function for the meson M. T^^^x) and T/^(^, x, u) are hard-scattering 
amplitudes, which are perturbatively calculable. The second term in Eq. (2) 
represents spectator contributions. 

This method works well for the case with two light mesons like vrvr or tiK [2,3], 
in which the final-state mesons carry large momenta. Interestingly enough, 
when there is a heavy quark in the final state such as 5 — i> D^n^ , this method 
still works when a spectator quark of the B meson is absorbed by, say, a D 
meson [2,4]. In Ref. [4], they calculated the ratio T(B -^ D+n-)/T(B -^ 
D+*7r-) or T(B -^ D+^*^p-)/T(B -^ D+'^*h-) since the Wilson coefficients 
were known at leading order at the time. The absolute branching ratios for 
B -^ D^^*'n^ were calculated recently since the Wilson coefficients are known 
at next-to-leading order accuracy [2,5].0 However, when the spectator quark 
is absorbed by a light quark, say, in i? — > Dtt'^, nonfactorizable contributions 
are infrared divergent, and the factorization breaks down. 

When we consider the decay B -^ J/ipK, it looks ambiguous at first sight 
whether we can apply the same method used in B —>■ tttt, or B -^ D^tt~, 
since the spectator quark in the B meson goes into a light K meson. However, 
what is special about J/ip is that the size of the charmonium is so small 
(~ l/asiTL^) that the charmonium has a negligible overlap with the {B, K) 
system, hence enabling the same improved factorization method in the decay 
B —^ J/ipK. As the explicit calculation shows below, the nonfactorizable 



In Ref. [2], Beneke et al. correctly pointed out the typos in Ref. [5]. 



contribution is infrared safe and the spectator contribution is suppressed in 
the heavy quark hmit. These facts indirectly justify the use of the improved 
factorization formula Eq. (2) in B ^ J/ipK. 

When the mass of the J/ip meson is not negligible, the light-cone wave function 
of the J/ip meson should include higher-twist contributions. The light-cone 
wave functions are obtained in powers of m^/E or Aqcb/E where E (~ rrib) 
is the energy of the J/ip meson. For B decays into two light mesons, the higher- 
twist contributions are negligible since they are of order Aqcb/E. However, for 
B -^ J/ipK, higher-twist contributions are important. Therefore we expect 
that the decay rate using only the leading, asymptotic wave function of J/ip 
will be smaller than the experimental result. 

Due to the nonzero mass of J/ip, we can think of several approximations with 
different limits. For example, we can take the infinite mass limit of the b quark 
in which mj, goes to infinity while m^ is fixed {m^/mi, ^0). In this case the 
result simply reduces to the case oi B ^ mr . But this limit is hardly realized 
in nature and it is not reasonable to compare the theoretical prediction in this 
limit with experimental data. We can consider another limit in which rub goes 
to infinity, while m^/mi, is held fixed. In this case, if m^ is heavy enough, 
it would be a better idea to employ a nonrelativistic wave function for J/ip. 
However, we can still use the light-cone wave function for J/ip. But we have 
to include the effects of higher-twist wave functions in order to compare with 
experimental data reasonably since these effects can be appreciable. 

From now on, we assume that, in the limit in which mf, goes to infinity, m^ is 
heavy enough to regard the size of the J/ip meson as small, but light enough 
to employ the leading-twist light-cone wave function for J/ip. It is difficult 
to satisfy this interesting, but hypothetical limit in reality, but we will not 
go further into the detail about this point. In this paper, we present the 
decay amplitude at next-to-leading order in the strong interaction, neglecting 
0(Aqcd/'"^6) corrections, employing the same technique to B ^ J/ipK as in 
B -^ Tin. 

One technical point is that we will keep the mass m^ of the J/ip meson in 
the hard scattering amplitude. In order for this to be theoretically consistent, 
we also have to consider the higher-twist wave functions for J/ip since higher- 
twist wave functions contain corrections of order m^/mb- However, the main 
issue is whether the hard-scattering amplitudes in the heavy quark limit are 
infrared finite at each order in Ug and l/nib- The effects of the higher-twist 
wave functions correspond to the expansion in l/rUb, and each wave function 
is, in general, an independent function. Therefore we have to verify that the 
hard-scattering amplitude convoluted with each independent wave function is 
finite. What we pursue in this paper is whether the hard-scattering amplitude 
convoluted with the leading-twist wave function for J/ip can be reliably cal- 



culated in the heavy quark hmit. The inclusion of higher- twist wave functions 
to compare with experimental data will be considered elsewhere. 



The effective Hamiltonian Hgg for B -^ J/ipK is written as 
^K^^^^-(CiOi + C2O2) - VtbV* 

i=3 



H,s = % (k,V;:(CiOi + C2O2) - VtbV* J2 c^o^) , (3) 

V ^ i=3 



where Cj are the Wilson coefficients at next-to-leading order evaluated at the 
renormalization scale /x. The relevant operators in Hgg are given by 

Oi = Sa'j''{l -75)Ca -0^7^(1 -75)^/3, 
02 = Sa'j''il - 75)C/3 ■ C^7m(1 " 75)&a, 

03 = s„7''(l - 75)6^ ■ ^ g^7^(l - 75)g/3, 

q 

04 = s«7''(l - 75)^/3 ■ E^/37^(l - 75)g„, 

q 

Oe = s«7''(l - l5)b(, ■ E^/37^(l + 75)g«, (4) 

q 

where a, (3 are the color indices and the sum over q runs over -u, rf, s, c and h. 

In calculating the decay amplitude, we introduce the vector and the tensor 
decay constants as [6] 

(J/V^|c(0)7^c(0)|0) = -2^m^e;, (5) 

(J/^|c(0)7^7.c(0)|0) = 2/J(e;p, - e>^), (6) 

where /^ is the decay constant which can be determined from the leptonic 
decay of J /il'- fl is the tensor decay constant arising from the tensor current, 
which formally depends on the renormalization scale. 

The distribution amplitudes on the light cone at leading-twist accuracy are 
expressed compactly as [7,8] 
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u^^in^j^i^) + n{fi>)pM^)\- (7) 

In Eq. (7), x is the momentum fraction of a c quark inside the J /ip meson, and 
the asymptotic wave functions 0,/,(a;), (p'^{x) for the J/ip meson are symmetric 



Fig. 1. Feynman diagrams for nonfactorizable contribution at order oig in 
B -^ J /'4}K. The dots represents the operators Oi {i = 1, 4, 6). 

functions under x ^^ 1—x. The asymptotic form of the distribution amphtudes 
0^(a;) and 4>^{x) is the same, given as 0^(a;) = (f>^{x) = 6a;(l — x). In the 
numerical analysis, we also consider the wave function of the form 0^(a;) = 
0^(x) = 6{x — |), which appeals to the intuitive expectation of the wave 
function. 

If we neglect the radiative corrections and the power corrections in Aqco/'^^fe, 
the factorized amplitude at leading order is written as 






We do not include the effects of the electroweak penguin operators since they 
are numerically small. Here A^ is the number of colors, m^ is the mass of J /ip, 
and e^ is the polarization vector of J /ip. Fi{q^) [i = 0, 1) are the form factors 
for B ^ K, which are given as 



V,= {K{p')\s^^b\B{pB)) 



ms-rriK 



i2 ^2 



{PB + P'), - ^ , ^ Pm Fi(p') + ^ , '' p.Foip'), (9) 

pZ J pZ 

where p = Pb — p' is the momentum of J/ip with p^ = m'^. From now on, we 
will neglect the kaon mass for simplicity. 

As we can see easily in Eq. (8), this amplitude depends on the renormalization 
scale fi since the Wilson coefficients depend on /i while the matrix elements 
of the operators are replaced by decay constants and form factors which are 
independent of /i. Therefore this amplitude is unphysical. However, if we in- 
clude the as correction to the amplitudes, it turns out that the /i dependence 
of the Wilson coefficients is cancelled and the overall amplitude is insensitive 
to the renormalization scale. 

Nonfactorizable contributions at order a.^ come from the radiative corrections 
of the operators Oi, O4 and Oq, and the relevant Feynman diagrams are shown 
in Fig. 1. The radiative corrections with a fermion loop do not contribute due 



to the color structure. For each operator Oi, O4 and Oq, if we add all the 
diagrams in Fig. 1 and symmetrize the result with respect to a; <-> 1 — a:, the 
infrared divergence of each diagram cancels and the remaining amplitude is 
infrared finite. One thing to note is that there appear imaginary parts in the 
nonfactorizable contribution, which are due to the final-state interaction. The 
strong phase can be calculated in the QCD-improved factorization and it is 
important in exploring the CP violation in nonleptonic decays. 

The decay amplitudes, in general, contain two terms proportional to Fq and 
Fi, which come from the matrix elements for B ^ K. However, if we use the 
equation of motion, we can find a relation between these terms, so that we can 
write the decay amplitude which depends only on Fq, or Fi. In the framework 
of the QCD-improved factorization with light-cone distribution functions for 
the mesons, we assume that all the quarks inside a meson are on their mass 
shell. Therefore we can safely use the equation of motion. For example, if we 
contract p^ to the current S7^(l — 75)6, we get 

J9^S7^(1 - 75)6 = ^bO- - l5)b = mbs{l + 75)6. (10) 



Similarly, if we contract the kaon momentum p'^ to the same current, we 
obtain 



s/(l- 75)^ = 0. (11) 



Here we neglect the mass of the strange quark, and accordingly the kaon mass, 
for simplicity. When we combine Eqs. (9) and (11), we obtain the relation 
between the two form factors Fq and Fi as 



ttirx „ , n, mi 



p'-V = ps- p'[{l - -f)Fr{p') + ^Foip' 



From this relation, the ratio Fq/ Fi is given by 

Fo(p2) p2 



0. (12) 



Fi(p^) m% 



This relation was also observed in calculating heavy-to-light form factors in 
the large energy limit [9]. In parameterizing these form factors using the QCD 
sum rule, this ratio is also valid to first order in p^ [10]. Note that, as p^ -^ 0, 
this ratio becomes 1 as it should be to remove the pole at p^ = 0. 

When we calculate the decay amplitude, some of the terms have the form 



(m| - ml)sf{l - 75)6 - 2mBe* ■ Pbs{1 + 75)6 

p ■ p's^*{l - 75)6 - e* ■ Pbs:^{1 - l^)h 
: 2{p ■ p't*^ - e* ■ py )s7m(1 - 75)& = 20^^57^(1 - ^^)h. 



(14) 



In deriving Eq. (14), we use the equation of motion with pb = p + p' , e* ■ p = 
and neglect the hght quark masses and the hght meson masses. The vector a^ 
satisfies the relation 



a ■ p' = 0, a ■ p = a ■ Pb = — e* ■ p'm\. 



(15) 



Therefore we obtain the relation 



a-V 



■Psm^ 



p2 



F,{p') + ^Foip' 

pZ 



(16) 



as a result of Eq. (12). With this identity, in calculating the decay amplitude, 
those terms proportional to e* ■ ^^57(1 + 75)6 can be replaced by the terms 
proportional to s^*(l — 75)^ or vice versa. 

The full amplitude for B -^ J/ipK is written as 



Gf 



iM = -2if^m^e* ■ VB—r= [KbV;s«2 - V't6V;;(a3 + a^)\ Fi{m, 



^Jj 



(17) 



where the coefficients a^ {i = 2, 3, 5) in the NDR scheme are given as 



Ci asCp^ I 



121n^ + /, + /,,), 



as = C, + ^ + ^^Cj-18 + 12ln^ + f, + fu), 



N An N 
Cq as Cp 

^ ^ A^ All N 



CJ-6 + 12ln^ + fj + fii 



The function // is given by 



fi= dx (f)^{x) 


1 — X 



r3('l-2a;), _ _ ,^ ^ 2r'^(l-x) 



1 — X 



In X — Sin + 3 ln(l — r ) + 



1 — r^x 



r^x^nn(l — r^) — inj 



( l-x X \4,2 '-^ V''"^-" " 

+ -, ^TTTT ^ v xuir X H ^ 

V(i_r2a;)2 {x_r\\-x)f' {\ - r\\ 



x)y 





Fig. 2. Feynman diagrams for the spectator contribution. 

+ -j^^ / dx 05;(a;)4r2[(- ^ _ _\ ^^^^2 

j^m^ J ^ ^^1 — ^^[1 — X) i — r'^x^ 

ln(l — r^) — i-K 



1 — r^fl — x) 



(19) 



where r = m^/mB- The function // is infrared finite even when r is nonzero. 
We organize the result such that it is proportional to Fi(m^), and all the 
terms involving r go to zero as r ^ 0. In this form, we can see clearly that 
the coefficients Oj in the limit r — i> is the same as the coefficients for, say, 
B ^ nir [2]. Note that the contribution from the tensor current vanishes as 
r ^ 0. This is trivial since the polarization vector e^ is proportional to the 
momentum p^^ in the massless limit and the tensor current does not contribute 
at all in this limit. 

/// is obtained from the the spectator contribution to B —^ J/ipK, which is 
shown in Fig. 2. This corresponds to the second term in Eq. (2). When we 
symmetrize the amplitude with respect to a; ^-* 1 — a;, it turns out that there 
is no infrared divergent part and the result is written as 



iMn = -2tUm^e* ■ ^^^(^5^:^ - Vtf,V*{C, - C,] 

^""^ A^^^^^TT^ J ^^^ J ^"^^ J ^"^' (20) 

where 0^, (px are the light-cone wave functions for the B meson and the K 
meson respectively. The spectator contribution depends on the wave function 
(pB through the integral 

\i^^^'-^. (21) 



Since 0b (0 is appreciable only for ^ of order Aqcd/ rriB, \b is of order Aqcd- 
Therefore /// is given by 

/.. = '4tJ4^T^ fd^"^ fdx't^ ]dJj^. (22) 



A^ Fi{rnV)rii?^ 1 — r^J S, I ^ I "" 



00 



The proof that the Q CD-improved factorization method works for B -^ J/ipK 
is sophisticated. There are two independent parts which depend on f^(j)^{x) 
and f^(f)j,{x). And each hard-scattering amphtude amphtude proportional to 
/^ or f^ is infrared finite. In both hmits of r -^ and r 7^ 0, the decay 
amphtudes are infrared finite, justifying that the method of the QCD-improved 
factorization can be apphed in B ^ J/ipK. This is the exphcit proof of the 
QCD-improved factorization for B -^ J/ipK at leading-twist accuracy, which 
was briefly mentioned in Ref. [2]. 

One important ingredient in the above argument is that it holds true only at 
leading-twist order. It means that we use only the leading-twist wave functions, 
and higher-twist wave functions are assumed to be suppressed. If E is the 
energy of the J/ip meson, it corresponds to neglecting those terms suppressed 
by {m^,AQCB)/E ~ (m^, Aqcd)/''^6- In reality, we expect that higher-twist 
effects are not negligible for appreciable values of r. When we consider higher- 
twist effects, we have to include the terms of order Aqcd/""^;) in the higher-twist 
wave functions as well as in the hard-scattering kernel. In order for the QCD- 
improved factorization to work, the higher-twist effects should also be infrared 
finite. 

For numerical analysis, we use the following input parameters: 



nib = 4.8 GeV, rric = 1.5 GeV, 
rriB = 5.28 GeV, m^ = 3.1 GeV, 
/^ =405 MeV, /b = 180 MeV, fx = 160 MeV. (23) 

We also choose A^^^ = 225 MeV, and A^ ~ 300 MeV. For the form factors in 
B ^ K, we use the values taken from Ref. [10]. 

Fi(mJ) = 0.606, Fo(mJ) = 0.418. (24) 



And the CKM matrix elements are expressed in terms of the Wolfenstein 
parameters [11] with A = 0.81 ± 0.06 and A = sin^c = 0.2205 ± 0.0018, 
and we fix them to their central values. For the wave function of J/ip, we 
use employ two kinds of the wave functions. One is the asymptotic function 
(pip{x) = 4>^{x) = 6a;(l — x), and the other is 0^(x) = (p'^{x) = 6{x — 1/2), 
which is more intuitive. 

In order to extract the magnitude of the tensor decay constant /J, we contract 
Eq. (6) with p'^, and use Eq. (5) and the equation of motion. Then we get the 
relation 

2mc(/x)/^ = m^/J(/x), (25) 



Table 1 

The coefRcients a^' at /j, = rrn, and mb/2 with different wave functions of J/V'- The 
values in the second and the third columns are the values with 0^(x) = 6x(l — x) 
and the fourth and the fifth columns are the values with the delta function. 



ai 


H = mb 


fi = mb/2 


/i = mh, delta 


^ 


= mh/2, delta 


02 X 10^ 


72.4-60.0i 


48.3-79.6i 


77.3 -GO.li 




54.9-79.8i 


03 X 10^ 


4.01+1.39i 


5.37+2.45i 


3.90 +1.39i 




5.17+2.46? 


05 X 10^ 


-4.95-1. 7H 


-6.03-3.97i 


-4.80-1. 71i 




-5.70-3.98i 



from which, we use 



'2mlU 



f^m, = ^^. (26) 



The numerical result of the coefficients Oj is summarized in Table 1. We can 
see that the coefficients Oj do not depend sensitively on the choice of the 
wave functions since the numerical values with different choices of the wave 
functions for J/ip are not much different. 

With these coefficients Oj, we can calculate the branching ratio. The branching 
ratio from experiment is given by [12] 

Bt(B -^ J/i)K) = (8.9 ± 1.2) X 10-^ (27) 

Our result shows that, for 4>ip(x) = 6a;(l —x), the branching ratio is 1.03 x 10~^ 
for fi = mi) and 1.05 x 10~^ for /i = 'mb/2, which are about eight times smaller 
than the experimental result. Though the theoretical treatment of the QCD- 
improved factorization in B -^ J/%pK is improved, the theoretical result still 
does not saturate the experimental result. This is expected since we do not 
include higher-twist effects which are supposed to be rather appreciable. In 
order to have a reasonable comparison, we need to include higher-twist effects, 
which is going to be pursued elsewhere. 

In Fig. 3, the dependence of the branching ratio on the renormalization scale 
is shown. The horizontal short-dashed lines in Fig. 3 show the experimen- 
tal branching ratio at la level. We also show the factorized result without 
nonfactorizable contribution. Of course, this quantity is unphysical and very 
sensitive to the renormalization scale, as we can see clearly in Fig. 3. On the 
other hand, for the full branching ratio calculated at next-to-leading order, the 
dependence of the branching ratio on the renormalization scale is very mild. 

Recent observation of CP asymmetry in B ^ J/ipKs [16] attracted some 
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Fig. 3. Dependence of the branching ratio on the renormahzation scale fi. The 
sohd hne uses the wave function 4>ii,{x) = 6x(l — x), the long-dashed line uses 
(t>^{x) = 5{x — i). The dotted line represents the result without the ag correction. 
The band with short dashes represents the experimental data at la level. 

interests in the possibility of new physics [17]. In this scheme in the standard 
model, there is no contribution to CP asymmetry in the decay amplitude 
since the CKM matrix elements involved here are all real. The CP asymmetry 
totally comes from B — B mixing, and we have not considered this issue here. 

In conclusion, we have considered the effect of nonfactorizable contributions in 
B -^ J /ipK in the heavy quark limit. In this limit, nonfactorizable contribu- 
tion can be systematically calculated from first principles using perturbation 
theory. It is shown that the decay amplitude is factorizable at next-to-leading 
order in the strong interaction and at leading order in Aqcd/""^?)- The nonfac- 
torizable contribution and the spectator contribution are infrared finite. 

In general, when the B meson decays into two mesons in which one is heavy, if 
the spectator quark of the B meson goes into the light meson, the decay am- 
plitude is not factorizable. The overlap of the heavy meson with the remainder 
is sizable and the soft gluon exchange gives a significant effect. However, when 
the heavy meson is the charmonium, the factorization method still works at 
leading-twist order. This physically means that the size of the charmonium 
is so small that the overlap with other mesons is small. Therefore the QCD- 
improved factorization method holds in B ^ J /il)K, in contrast to naive 
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expectations. 

In spite of the theoretical improvement mentioned so far, the theoretical 
branching ratio is about eight times smaller than the experimental result. 
This is typical in class II decays in which there is color suppression in the 
lowest order contribution. And the as corrections give a significant change in 
magnitude and a significant strong phase compared to the leading-order result 
corresponding to naive factorization. 

Another reason why the theoretical result is small compared to the experimen- 
tal data is that we used only the leading-twist wave functions. In a realistic 
case, higher-twist effects are not negligible. Khodjamirian and Riickl [13] con- 
sidered nonperturbative effects in B —^ J /ipK using the QCD sum rules, 
and they found that nonperturbative effects including the higher-twist wave 
functions could indeed be large. Also in the inclusive production of J /ip., the 
naive factorization does not explain the large production rate. It was suggested 
that the octet contribution would enhance the rate [14,15]. This contribution 
comes from the nonleading Fock state \gcc) state, which constitutes higher- 
twist effects, and it may not be negligible in the decay B -^ J /ipK either. 
However, the method of the QCD-improved factorization can be safely ap- 
plied to -B — s> J /ipK at leading-twist order and it is remarkable that we obtain 
a very significant correction to naive factorization. 

Note added: While this paper has been written, Ref. [18] appeared, in which 
the decay amplitude for B -^ J /ipK are computed including chirally enhanced 
contributions. Their result looks different from ours, but the definition of // 
and gi are rather ambiguous since they can be switched back and forth us- 
ing the equation of motion discussed in this paper. Furthermore, they put 
ferric/ f^tpm^ = 2x^ using the equation of motion. That is why there is no /J 
in their result. However, when all is taken into account, both results are the 
same. 
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